Mixed-state microwave response in superconducting cuprates. 
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We report measurements of the magnetic-field induced microwave complex resistivity in 
REBa2Cu307_5 thin films, with RE = Y , Sm. Measurements are performed at 48 GHz by means 
of a resonant cavity in the end-wall-replacement configuration. The magnetic field dependence is 
investigated by applying a moderate (0.8 T) magnetic field along the c-axis. The measured vortex 
state complex resistivity in YBa2Cu307_5 and SmBa2Cu307_5 is analyzed within the well-known 
models for vortex dynamics. It is shown that attributing the observed response to vortex motion 
alone leads to inconsistencies in the as-determined vortex parameters (such as the vortex viscosity 
and the pinning constant). By contrast, attributing the entire response to field-induced pair break- 
ing leads to a nearly quantitative description of the data. 
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I. INTRODUCTION 



The microwave response is a source of important infor- 
mation in superconductors. In high- Tc superconductors 
(HTCS) various microwave techniques have been used to 
get information, among the others, on the symmetry of 
the order parameter, on the vortex parameters such as 
the vortex viscosity and pinning frequency and on the 
temperature dependence of the superfluid fraction (via 
the measurement of the temperature dependence of the 
London penetration depth). In short, microwave mea- 
surements are a very powerful tool to study the super- 
conducting state. 

While YBazCuaOy-a (YBCO) and Bi2Sr2CaCu208+:r 
(BSCCO) have been the subject of intensive experimen- 
tal investigation at microwave frequencies, other HTCS 
did not receive the same attention. In particular, a very 
few reports dealt with the vortex-state microwave re- 
sponse in rare-earth substituted 123 compounds, mainly 
on Gd-substituted materials^. In fact, RE-substituted 
123 compounds present interesting features in the vortex 
state: in particular, they often give enhanced irreversibil- 
ity lines and in some cases a slightly higher Tc- Due to 
their very similarity to YBCO, they are not thought to 
present novel features in the basic mechanisms govern- 
ing the superconducting properties. Aim of this paper is 
to present a compared study of the microwave response 
of YBCO and SmBazCugOy-i (SmBCO) in the vortex 
state in moderate fields. It will be shown that some quan- 
titative and qualitative differences exist between the field 
dependence of the complex resistivity in the two materi- 



als. However, the analysis of the data shows that those 
differences cannot be related only to changes in the vor- 
tex dynamics (e.g., different pinning). Indeed, a close 
analysis of the data shows that vortex motion cannot be 
the main contribution to the observed resistivity in any of 
the two materials. Rather, we find that the resistivity of 
both YBCO and SmBCO can be qualitatively described 
assuming that the response is mainly determined by mag- 
netic field induced pair breaking. Within this frame, the 
differences between the two materials can be ascribed to 
different quasiparticle scattering time below Tc. 
We recall the main phenomena lying at the base of the 
microwave response in the vortex state. In principle, the 
response is dictated by vortex motion, described by the 
vortex motion complex resistivity Pvm, and by charge 
carriers complex conductivity, a. The latter originates 
from superfluid as well as normal carriers. 
A general frame that includes those contributions is the 
CofFey-Clem theory^. It is noteworthy that the theory, 
while developed for a particular case of vortex motion 
(vortex in a sinusoidal potential submitted to the Lorentz 
force and to a stochastic force), is not constrained to 
a specific potential shape seen by the vortices. Given 
the vortex motion resistivity pmn and the charge carriers 
complex conductivity a — ui — ia2, the resulting expres- 
sion for the total complex resistivity can be written as: 



(1) 
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When the field dependence of the pair breaking (affect- 
ing a) is negligible, and for T sufficiently below such 
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that (72 ^ CTi, one has for the field induced change of the 
resistivity: 



Ap = p{H, T) - p(0, T) = p,miH, T) 



(2) 



since the vortex motion resistivity is zero by definition in 
zero field. By contrast, when the vortex motion resistiv- 
ity is vanishingly small, equation 1 reduces to: 



1 



a{H,T) 



(3) 



The vortex motion approximation alone has been often 
used for the analysis of the microwave response in the 
vortex state, ^ sometimes including a field-independent 
contribution from the superfluid.^"^ In this case, from 
the data one can estimate vortex parameters such as the 
vortex viscosity, the depinning frequency and the pin- 
ning constant and the activation energy. By contrast, 
it has been shown that in BSCCO films the imaginary 
conductivity could be well described'' by expression 3. 
In this latter case characteristic parameters determining 
the response are the quasiparticle scattering time and the 
pair-breaking field. 



II. EXPERIMENTAL SECTION 

The SmBa2Cu307 film was grown on 1 cm x 1 cm 
LaAlOs substrate by planar high oxygen pressure dc 
sputtering technique. The roughness of ^ 20 nm on 1 /im 
X 1 /zm area was measured by AFM. The YBa2Cu307 
film was grown by inverted cylindrical magnetron sput- 
tering on 1 cm X 1 cm Ce02 bufi^ered YSZ substrates. 
In both samples, out-of-plane c-axis orientation was as- 
sessed by rocking curves, that showed FWHM w 0.2°. 
Film thickness was d ~200 nm. $-scan showed good epi- 
taxiality. YBCO and SmBCO showed ~ 88 K and 87 
K, respectively. Further details on sample preparation 
have been reported elsewhere.*'^ 

The microwave response was measured by the end-wall 
cavity technique at 48.2 GHz in zero dc magnetic field, 
and with an applied field (along the c-axis) /zoi? <0.8 T. 
The cylindrical cavity resonated in the TEqh mode.^° 
In this configuration the microwave field probed a cir- 
cular area, corresponding to the diameter of the cavity, 
^ 8 mm, centered on the center of the film. We note 
that in this configuration the edges of the film are not 
exposed to the microwave field, so that our measure- 
ments are not affected by detrimentals edge effects (un- 
like, e.g., dc and microwave measurements on striplines). 
Due to the field distribution of the excited mode, the 
peak microwave field /io-ff^uj ^ 10 /iT is reached in an 
annular area of 4 mm diameter centered on the center of 
the film. Thus, for all practical purposes, the amplitude 
of the microwave field _ff^^ ^ H, so all the field effect 
should be regarded as coming from the dc applied field. 
Field-induced variations of the quality factor and 



of the resonance frequency Av were measured at several 
temperatures in the range from 67 (69) to 90 K, result- 
ing in reduced temperature ranges above T/Tc = 0.76 
(0.79) in YBCO (SmBCO). Data have been collected ei- 
ther in zero-field-cooled, field-cooled, and on direct and 
reverse field-sweeping. Wc observed no hysteresis, with 
the possible exception of the data below ~70 K, wh(Te 
a small hysteresis, smaller than 10% of the overall sig- 
nal is occasionally present on reverse field sweeps. In 
this case, data for field sweeps after ZFC are considered. 
We note that rather high operating temperatures and 
the strong demagnetization present in our arrangement 
(the dc field is perpendicular to the square film) lead to 
a very small value of the first penetration field, and in 
particular to B ~ iiqH even at low fields. However, 
for conservative reasons, we will not focus the discussion 
on the low field data (below 10 mT). Figure 1 reports 
the temperature variations of the cavity parameters re- 
flecting the superconducting transition in zero magnetic 
field in YBCO (left panel) and SmBCO (right panel). As 
can be seen, a substantial background is present in the 
measurement of the frequency shift, indicating relevant 
mechanical relaxations of the metallic resonator. Even 
if the background was rather reproducible, we have cho- 
sen to work with field variations at fixed temperatures. 
In this case, the variation of the effective impedance 
does not depend on the background (but only on a ge- 
ometrical factor G, which acts scale factor), and 
and Ap directly yield the field induced changes of 
the effective surface impedance of the thin film, accord- 



ing to: G \^A^-i2^^ = AZffiH) ^ [ApiH)]/d = 

[Api{H) + iA/92(H)] /d, where p is the complex resistiv- 
ity. Here, we have made use of the thin film approxima- 
tion, extensively examined in Ref. 12. 
In figure 2 we show the effect of the application of a 0.7 T 
magnetic field in YBCO and SmBCO, respectively. The 
field-induced dissipation has a similar behaviour in both 
materials: [/9i(T, 0.7T) — /3i(T, 0)] /d increases with in- 
creasing temperature, and then vanishes at the transition 
temperature, as qualitatively expected. The variation 
of the imaginary part [p2{T, 0.7T) - p2(T, 0)] /d exhibits 
a quantitatively different behaviour in the two materi- 
als: in YBCO it is positive but small, apart the region 
very close to the transition where, at a temperature Tq, 
it becomes negative and eventually approaches zero for 
T Tc- In SmBCO the qualitative behaviour is similar, 
with a larger low temperature value and a much lower Tq . 
A deeper investigation is made by measuring the field- 
sweeps at fixed temperature. In figures 3 and 4 we show 
typical field sweeps, Ap{H)/d = [p{H,T) - p{0,T)] /d, 
in YBCO and SmBCO, respectively, at temperatures 
chosen in both cases a few Kelvins below Tq. The data 
here reported are representative of the behavior observed 
in other SmBCO and YBaCuO films measured under 
the same conditions and in similar temperature and field 
ranges. It is seen that the field dependences of the mi- 
crowave resistivity in YBCO and SmBCO are slightly 
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different. In particular, in the present field range the 
field variation of both Api and A/92 in SmBCO presents 
a more pronounced downward curvature than in YBCO. 
We discuss in the following the possibile reasons for the 
similarities and differences between YBCO and SmBCO. 



III. DISCUSSION 

When wc come to disciiss the data, we arc faced with 
the increase of the imaginary part in even moderate fields 
at low enough temperatures (see figures 2,3,4). Since the 
role of pinning reveals itself mostly on the imaginary part, 
one might be tempted to assign the quantitative differ- 
ence between YBCO and SmBCO to a different pinning, 
stronger in SmBCO. Following this quite conventional 
view, we first discuss the data in terms of simple vortex 
dynamics, assuming that the observed resistivity is en- 
tirely due to vortex motion (VM). To this aim, we take 
the simplest model for as early proposed by Gittle- 
man and Rosenblum^^ (GR): 



^qB 1 + i^ 



(4) 



where r] is the unit length vortex viscosity, u>p is the 
depinning angular frequency and kp = ajp-ri is the unit 
length pinning constant. This analysis has been widely 
used in various HTCS, see the review in Ref. 3. This 
model, based on the motion of independent vortices, is 
expected to be more and more accurate with increasing 
frequency: at high operating frequencies vortex displace- 
ments induced by the microwave field are much smaller 
than the intervortex mean spacing, and complications 
arising from the collective nature of the vortex matter 
are less and less relevant. In fact, it has been exper- 
imentally shown^'' that in YBCO a dynamic crossover 
from a glassy behaviour to the independent vortex mo- 
tion takes place above a few GHz. We note that the GR 
model explicitly predicts that the imaginary resistivity 
should increase with the application of the field. This 
result is common to other VM based models, based on 
less restrictive assumptions. It is then apparent that, 
at least in a temperature range of several Kelvins below 
Tc (above Tq), the data cannot be explained on the ba- 
sis of VM alone. At low enough temperatures, applying 
in a straightforward way the GR model one would get 
directly from the data the depinning frequency ujp/2TT, 
the vortex viscosity rj and the pinning constant kp. In 
figures 3 and 4 (right panels) we report the calculated 
vortex parameters in YBCO and SmBCO. We immedi- 
ately note that in both materials the vortex viscosity and 
pinning constant increase with the field, while the depin- 
ning frequency increases in YBCO and is approximately 
constant in SmBCO. The field dependence of the so- 
obtained r] and hp is not easily explained. In particular, 
one would expect a nearly constant vortex viscosity, as 



experimentally determined by multifrequency measure- 
ments at lower temperatures® and by high-field swept- 
frequency Corbino disk measurements,^^ and a weak field 
dependence of kp. Exotic field dependences of the vortex 
viscosity, increasing with the field as ^ Vb, have been 
observed in Bi2Sr2Cu06 and tentatively explained^® in 
terms of a different field dependence of quasiparticle re- 
laxation time in a d-wave superconductor, appearing at 
high fields when the intervortex distance becomes smaller 
than the mean free path. This picture does not apply 
easily to our case, since the field range is very differ- 
ent. In addition, this interpretation would not give an 
explanation of the field increase of kp (we remark that 
the anomalies reported in Ref. 16 were confined to r]: 
kp was found constant at low fields, and decreasing at 
higher fields). We conclude that this explanation is at 
least questionable at low fields. 

Granularity is sometimes indicated as a possible dom- 
inant source for the losses in the microwave response 
in superconducting films. Manifestations of granular- 
ity include weak-links dephasing,^^ Josephson fluxon 
(JF) dynamics^* and, as recently studied, Abrikosov- 
Josephson fluxon (AJF) dynamics.^'' Weak- links dephas- 
ing is charaterized by a very sharp increase of the dissi- 
pation at dc fields of order or less than 20 mT, accom- 
panied by a strong (and sometimes exceptionally strong) 
hysteresis with increasing or decreasing field. However, 
those effects are relevant in large-angle grain boundaries, 
such as those found in pellets and granular samples, and 
are not observed in good thin films. In fact, we did not 
observe none of the above mentioned effects in our films, 
so we conclude that weak-links dephasing does not affect 
our measurements. Josephson fluxon dynamics has been 
studied essentially in relation to nonlinear effects, due to 
the short JF nucleation time. If however one assumes 
that a dc field has the same effect as a microwave field, 
the qualitative properties of JF dynamics would be barely 
distinguishable from conventional Abrikosov fiuxons dy- 
namics, since an equation like Eq.4 would hold,^* with 
the noticeable difference of a larger pinning frequency 
(due to the small viscosity, due in its turn to the insu- 
lating core). Reactive response should then dominate, 
while in our measurements (see figs. 2,3,4) the dissipa- 
tion (Api) is always larger than reactance (A/92)- We 
are led to conclude that JF dynamics, even if possibly 
present, is not the main mechanism driving the measured 
microwave response. 

Another possibility is the role played by Abrikosov- 
Josephson fiuxons. Such flux structures nucleate 
along small-angle grain boundaries, as those presumably 
present in our oriented films. The ac response is expected 
to saturate at fields larger than a characteristic field 
Hq, whose estimate spans orders of magnitude, ""^^ in the 
range 0.1-10 T. In YBCO bicrystals, dc measurements^^ 
yielded iJo < 2 T above T = 70 K. However, it is 
noteworthy that in defective thin films of other cuprates 
(Tl2Ba2CaCu208+a:) the 8.5 GHz dissipation ascribed to 
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AJ fluxon motion is shown to saturate at (microwave) 
fields as low as ij-oHq ^ 1 mT,^^ while in better films this 
kind of response is not present at all. Nevertheless, allow- 
ing for a large Hq, the explicit expression for A J fluxon 
motion complex resistivity^^ yields an initial magnetic- 
field increase in the resistivity as ^ \f^y whereas the 
reactance initially decreases with the magnetic field as 
- + c{T) with - i < c < i . It appears that 

our data, where the reactance increases or decreases with 
the field dependending on the temperature, are not easily 

described by this only mechanism. 

Summarizing, none of the mechanisms shortly recalled 
above appears to describe the main qualitative features 
of our data, namely: (i) the sublincar variation of the 
real and imaginary parts of the resistivity, ( ii) the change 
of sign of the variation of the reactance with increasing 
temperature. It is not excluded that a properly tailored 
combination of several of the above mechanisms could 
describe our data. Nevertheless, in order to look for a 
single possibly dominant mechanism in the microwave 
complex response to a dc magnetic field we now com- 
ment on a possible alternative explanation for our data. 
In fact there is quite robust evidence, in other HTCS, 
that in the vortex state the quasiparticle increase has a 
relevant role, at least in the imaginary response.'' This 
phenomenon is well explained in terms of a semiclassi- 
cal theory of a superconductor with line of nodes in the 
gap^^. In particular, the decrease of superfluid fraction is 
predicted to follow a B" law, with a = ^ in clean d-wave 
superconductors. In this case the charge carrier conduc- 
tivity can be described by the two-fluid model, and cast 
in the form: 

a{H,T) = aRQXniT,B) - ia-/oXs(r, B) (5) 

where a;„(T, B) and Xs{T, B) = 1 — x„(T, B) are the tem- 
perature and field dependent normal and superfluid frac- 
tional densities, and Xs{T,B) = .t^(T,0) [1 - [B / Bpb)"] 
where iJpf, is a pair breaking field. It has been shown^'^ 
that in clean superconductors Bph oc Bc2{T), but devia- 
tions from this behaviour are possible. 
As a very crude approximation, we can assume as a lim- 
iting case the vortex motion does not give a significant 
contribution to the resistivity. Taking into account that 
{B/Bpb)°' <C 1, one can show-^** that in this case the field 
variation of the complex resistivity can be written as: 

Ap(//,r) = [ai(T;s)+ia2(T;s)] (^-^^ (6) 

where the dependence on the single parameter s = 
has been made explicit. By considering the tempera- 
ture dependent data in figure 2, we do not need to as- 
sume any specific field dependence. However, we remark 
that in SmBCO the data are well described with a = ^ 
(Ref.^^). Due to the oversimplification of the model, 
for the temperature dependences involved in Equation 



6 we take conservatively the most simple: the super- 
fluid fraction Xe{T,Q) = (l-i^), and Bpfe oc Bc2{T), 
so that Bpb = Bpbo (l - f^), with t = T/Tc. In order to 
gain qualitative information, we do not attempt to in- 
sert some temperature dependence in the parameter s. 
We stress that, with this choice, s is not a free parame- 
ter: its value is determined by the requirement that the 
calculated curve of Ap2 changes sign at the experimental 
temperature Tq. The only free parameter is the scale fac- 
tor given by the pair breaking field. The model is surely 
oversimplified, but has rather strong constraints: in par- 
ticular the shapes of the two calculated curves Api{T) 
and Ap2{T) are determined, and only a common scale 
factor can be adjusted. The resulting calculated curves 
for the experimentally applied field Ba = 0.7 T are re- 
ported in figure 2. As can be seen, the overall shape of the 
experimental curves arc reproduced. Prom the scale fac- 
tors we estimate ( ^ 0.05 in YBCO and SmBCO. 

One also gets s ~ 0.1 and 0.2 in YBCO and SmBCO, 
respectively. We note that by applying the conventional 
two-fluid model, a = ^ [torXn — ixs], one has s = ljt 
where r is the quasiparticle scattering time. Taking into 
account the measuring frequency, one has r ~ 0.35 ps 
and 0.7 ps in YBCO and SmBCO, respectively, that com- 
pare well to, e.g., 0.2 ps at «80 K as obtained from mi- 
crowave measurements in YBCO crystals^^ and to 0.5 ps 
at s=;80 K as obtained from millimeter-wave interferome- 
try in YBCO film.^^ We note that the model appears to 
describe better the data on SmBCO than on YBCO: this 
is in fact consistent with the much stronger curvature of 
the field-sweeps data in SmBCO, since the proposed field- 
induced pair breaking mechanism is naturally sublinear 
in superconductors with nodes in the gap. In fact, in 
SmBCO the small difference between the calculated pair 
breaking response and the experimental data can be well 
described by vortex motion in the free-flow limit. ^'''^^ In 
this respect, the discrepancy between the calculated pair- 
breaking and the experimental data, as well as deviations 
from if 2 dependence of the field sweeps, is most likely 
due to motion of flux lines, either JF, AJF or Abrikosov 
fluxons. We believe that the successful application of the 
oversimplified pair-breaking model in the qualitative de- 
scription of the vortex-state resistivity in RE(Y,Sm)BCO 
points toward a substantial contribution of pair breaking 
by the magnetic field, in agreement with the existence of 
lines of nodes in the gap, while additional mechanisms, 
such as fluxon motion, are most probably also present. 

IV. CONCLUSION 

We have presented data for the microwave resistivity 
in REBa2Cu307_5 thin films as a function of the tem- 
perature and magnetic field. The field-dependence of the 
complex resistivity exhibits quantitatively different be- 
haviours in the two materials. The analysis in terms of 
vortex motion alone, in particular ascribing the quan- 
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titativc differences to a different pinning, would imply 
very exotic field dependences of the vortex parameters, 
so that it appears that alternative interpretations have 
to be found, including a possible enhanced role of the 
field-induced pair breaking. By ascribing, as a limiting 
case, the entire response to the field-induced pair-braking 
we obtain semi-quantitative fits of the data at 0.7 T as 
a function of the temperature, with a largo cjuasiparti- 
cle scattering time below Tc in agreement with estimates 
given by several groups. Additional measurements might 
be required in order to fully understand the interplay 
between fluxon motion and pair breaking. 
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FIG. 1. Temperature dependence of the coeffcients of the 
resonant cavity incorporating as an end wall the supercon- 
ducting film. Main panels, quality factors; insets, resonant 
frequency. Left panel: YBCO film; right panel, SmBCO film. 




FIG. 2. Variation of the resistivity with the application of 
a 0.7 T magnetic field. Left panel reports data on YBCO, 
right panel in SmBCO. Full symbols; shift of the real part, 
[pi(T,0.7T) - pi(T,0)] /d. Open circles; shift in the imag- 
inary part, [^2(7^, 0.7T) — p2{T,0)] /d. Lines are calculated 
curves with Eq.6 with a common scale factor as the only free 
parameter. 



FIG. 3. Left panel; field induced change in the real (full 
circles) and imaginary (open circles) part of the complex 
resistivity. Right panels; vortex parameters as calculated 
from data in YBCO accordingly to the conventional Git- 
tlemam-Rosenblum model. From top to bottom; depinning 
frequency uJp/2t:, unit length vortex viscosity 7], unit length 
pinning constant kp. The strong field dependence of the 
so-calculated vortex viscosity and pinning constant cannot be 
easily justified. 
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FIG. 4. Same as in Figure 3, but for the SmBCO film. 
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